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£Nj ■ Abstract 

The dark energy issue is attracting the attention of an increasing number of physicists 
all over the world. Among the possible alternatives to explain what as been named the 
"Mystery of the Millennium" are the so-called Modified Theories of Gravity. A crucial test 
for such models is represented by the existence and (if this is the case) the properties of their 
}_( ' black hole solutions. Nowadays, to our knowledge, only two non-trivial, static, spherically 

symmetric, solutions with vanishing cosmological constant are known by Barrow & Clifton 
(2005) and Deser, Sarioglu & Tekin (2008). The aim of the paper is to discuss some features 
\ of such solutions, with emphasis on their thermodynamic properties such as entropy and 

■ temperature. 

00 ■ I Introduction 

, Since the discovery by Riess and Perlmutter and respective collaborators [T] E] that the universe 

is — against any previous belief — in an accelerating epoch, the dark energy issue has become the 
"Mystery of the Millennium" [3] . Today, dark energy is probably the most ambitious and tantalizing 
field of research because of its implications in fundamental physics. That the dark energy fluid 
has an equation of state index w very close to minus one represents an important point in favour 
of those who propose to explain dark energy in terms of a cosmological constant, A. Still, a non- 
vanishing cosmological constant does not exhaust the range of models that have been proposed 
so far in order to solve the aforementioned issue. This is justified, in part, by the whole sort of 
well-known problems raised by the existence of a strictly positive cosmological constant. 

On the other hand, it is well accepted the idea according to which general relativity is not the 
ultimate theory of gravity, but an extremely good approximation valid in the present day range 
of detection. It basically comes from this viewpoint the input to so-called modified theories of 
gravity which nowadays enjoy great popularity (c/. [H El El El El El EH] for a review). Without 
any claim for unification, such models propose to change the Einstein-Hilbert Lagrangian to a 
more general form able to reproduce the same general relativity tests on solar distance scales and 
further justify both inflationary and current acceleration of the universe. 

The original idea of introducing a correction to the Einstein-Hilbert action in the form of 
f(R) = R + R 2 was proposed long time ago by Starobinsky [11J in order to solve many of the 
problems left open by the so-called hot universe scenario. This, in turn, had the consequence of 
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introducing an accelerating expansion in the primordial universe, so that the Starobinsky model 
can be considered as the first inflationary models. The recent interest in models of modified gravity 
instead, grew up in cosmology with the appearance of |12[ [T51 [TIL 

The mathematical structure of /(i?)-theories of gravity and their physical properties (e.g., 
asymptotic flatness, renormalizability, unitarity) have been an exciting field of research over the 
last four decades; a small but significant trace of which is represented by |15IITIfl[T71[T81[T!Jll2T)ll21| . 

The arena of models is in principle infinite while departures from Einstein's theory are most of 
the times all but minimal. Of crucial interest is, of course, the existence and, if this is the case, the 
properties of black holes in modified gravities. It is quite easy to find the conditions allowing the 
existence of de Sitter-Schwarzschild black holes (see, for example [22J for f(R) modified gravity, 
[2"5] for Gauss-Bonnet modified gravity, and [241 [231 121)1 12"7] for related topics). 

Here, we are interested in non-trivial and static black holes solutions. However, the number 
of exact non-trivial static black hole solutions so far known in modified theories of gravity is 
extremely small: just two, both spherically symmetric. They have been obtained by Barrow & 
Clifton (2005) in a modified theory of the type f(R) = R 1+& with S a small real parameter; and 
by Deser, Sarioglu & Tekin (2008) by adding to Einstein-Hilbert Lagrangian a non-polynomial 
contribution of the type VC 2 , with C a bcd being the Weyl tensor. 

These black hole solutions are not expected to share the same laws of their Einsteinian coun- 
terparts: for this reason, following [3 1J . we shall refer to them as dirty black holes. Some of the 
physical quantities one would like to address to dirty black holes are their mass, the horizon en- 
tropy, their temperature and so on. Thanks to the large amount of work carried over in the last 
decade, we can firmly say that the issue of entropy and temperature of dirty black holes represents 
a well posed problem [30] ; a nice and recent review on the entropy issue associated with f(R) 
gravity models is [32J, where a complete list of references can be found. Here, we only mention 
[551 l2"2l |2"5] , However, with regard to the mass issue, all considerations still lay on a much more 
precise ground. 

In the present paper we shall work in units of c = G = h = Ub = 1< The organization is as 
follows: in §2 we review the Deser-Sarioglu- Tekin solution and compute entropy and temperature 
for such black hole; in §3 we do the same for the Clifton-Barrow solution. In the Conclusions 
we address the difficulties faced trying to define meaningfully the concept of mass for dirty black 
holes. 



II The Deser- Sarioglu- Tekin solution 

Let us start by recapitulating the Deser-Sarioglu- Tekin solution [55] • The authors start from the 
action 

I DST = / (fxJ^g (R + VSaVc 2 ) + Boundary Term (II.l) 

16tt J M ■ \ J 

where a is a real parameter and C 2 := C ab cd C cd ab is the trace of the Weyl tensor squared. Looking 
for static, spherically symmetric solutions of the type, 

dr 2 

ds 2 = -~a(r)b(r) 2 dt 2 + — + r 2 (d8 2 + sin 2 Odcf) {11.2) 
a{r) 



the action (III. Ill becomes 



I D ST[a{r)Xr)] = \ J dt J dr [(1 - a)(ra{r)b' (r) + b(r)) + 3aa(r)b(r)] . (II.3) 

Imposing the stationarity condition 5I[a(r),b(r)] — gives the equations of motion for the un- 
known functions a(r) and b(r). 

(l-a)rb'(r)+3ab(r) = 
(1 - o-)ra'(r) + (1 - 4a)a(r) = I- a. (II.4) 

According to a, the space of solutions of (|II.4|> can be different, in particular: 
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• cr = corresponds to Einstein-Hilbert action. In fact, a(r) = 1 — - and b(r) = k and for c, fc 
positive constants, the Schwarzschild solution of general relativity is recovered; 

• cr=l: only the trivial, physically unacceptable, solution a(r) = = b(r) exists; 

• cr = -t\ then, for some positive constants k and tq: 

a(r) = bx(-) and bir) = - ; (II.5) 
v r / r 

• In all other cases, the general solution to (|II.4|) turns out to be 

1 (7 14 f f \ 

a(r) = cr~ i-^ and 6(r) = ( — ) ° , (H-6) 

l — 4cr \fe/ 

for some positive constants c, k. 

The constants fc, and k appearing in b(r) are removable by time re-scaling. Notice also that, in 
(|II.5|) . goo and gn go to zero as r — > oo so that the model is unphysical. For this reason, we shall 
mainly concentrate on the solution (|II.6|) parametrized by all the a ^ 0, 1, \ . 
In order to treat pi.6jl . let us introduce the parameter p(a) := ^Er so that the metric becomes 

ds 2 = -(p - cr- 1 ,) (j) 2{ljE) dt 2 + d ^— r - + r 2 (d9 2 + sin 2 9dtf) . (11.7) 

\k/ (p — cr p) 

For p < 0, or ^ < a < 1, a(r) — — (\p\ + cr^) < for all r, that is, the parameter region 
j < a < 1 needs to be excluded to preserve the metric signature. As regard the asymptotic 
behaviour of (III. 71) . we see that: 

• for p>lor0<cr<j, we have that goo — > and gn — > - as r — > oo; 

• for < p < 1 or (T £ (— oo, 0) U (1, +oo), we have that goo oo and gn — > ^ as r — > oo. 

As noted by Deser et al. the fact that the asymptotics of goo and gn differ means that the 
equivalence principle is violated: something which is intimately related with the difficulty of 
defining a "mass" in this theory [28J . 

Looking at the solution (|II.7|) . we see that the hypersurface r = r# := f defined by the con- 
dition a(rn) — behaves as a Killing horizon with respect to the timelike Killing vector field £ a . To 
prove this, let us define a complex null tetrad j7 a , n a , m a , fh a } for the metric (|II.7|) according to the 
following 
rules [33] : 

1. l a is s.t. on the horizon 

l a H = £ a ; (11.8) 

2. The normalization conditions hold 

l-n = -l & m-m = l; (II.9) 

3. All the other scalar products vanishes. 

Since the metric (|II.7|) is not asymptotically flat, it is not clear at all what is the right normalization 
for £ a . Assuming £ a = \ d t a , A G K+, it's not difficult to see that 

r = (A, A a(r)b(r), 0,0) , 
1 1 



m a = 0,0 



2A a(r)b(r) 21 2Xb(r) 
i 1 



,0,0 



\/2r ' \[2r sin 
% 1 
V^r' ^rsin^ 



0,0,-4-. /s 1 . J - (n.io) 
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satisfy the list of conditions to form a complex null tetrad. As a consequence, for example, the 
metric can be re-written as g ab = —2l^ a n b ) + 2mr a fhb) . The null expansions are, by definition, 



e. 



V„n a + n a l b V a n b + n b l a V a n b 



0+ := VJ a + l a n b VJ b + l b n a VJ b 



2 A 



\rb(r) 
a(r)b(r) 



(11.11) 



Thus, in-going light rays always converge (0_ < for all r > 0); out-going light rays, instead, 
focus inside the horizon (0 + < as r < r#), diverge outside it (0 + > as r > rjj) and run in 
parallel at the horizon (0+|# = 0). When they are slightly perturbed in the in-direction (that 
is, along n), the out-going null ray is absorbed inside the horizon th as it is confirmed by the 
fact that the in-going Lie derivative C n Q + \H = — \- < is everywhere negative. Computing the 

convergence (g := —m a fh b V b l a ) and the shear (fi := —m a m b V b l a ) of the null congruences at the 
horizon we can immediately check they vanish, as expected for any Killing horizon. The Killing 
surface gravity 

k h := -l a n b VJ b \ H - \a'(r) H b(r) H , (11.12) 

turns out to depend by the normalization of the Killing vector £ a . In order to fix A, we may 
implement the conical singularity method. To this aim, let us start by the Euclidean metric 



ds% 



dr 2 
W{r) 



+ V(r)dr 2 + r 2 dfl 2 



where we suppose that both V(r) and W(r) have a structure like 

V(r) — (r — f)v(r) & W(r) = (r — f)w(r), 



(11.13) 



(11.14) 



with v(r), w(r) regular for r > f. f may be identified with some type of horizon close to which we 
are interested in the behaviour of the metric. 



O 2 , 



with C a constant we are going to fix very soon. 



ds% 



1 



w(r) 



dr' 



- + (r — r)v(r)w(r)dT 2 



+ (r + (x 2 ) 2 dfl 2 



4C 



w(f) 



'—dx 2 + Cv(r)x 2 dr 2 + f z dfl 



Let us choose ( = w(f)/A. the Euclidean metric takes the form 



ds 2 P rs dx 



l dn 2 



x < 1. 



(11.15) 



(11.16) 



(11.17) 



(|II.17P shows that close to the horizon (r ps f or x <€. 1) the metric factorizes into IC2 x § 2 f : IC2 
being the metric of flat two-dimensional metric on behalf of identifying x with the polar distance 
and r with the angular coordinate. However, IC2 is regular if and only if 



y/v(f)w(r) y/v(r)w{f) 
t *~ r + Zn 



or, in other words, 



An 



v{r)w(f) 



(11.18) 
(11.19) 



4 



P representing the (unique) r-period which allows to impose a smooth flat metric on IR 2 . 



In Quantum Field Theory, the KMS propagator exhibits a periodicity in time when the sys- 
tem is at finite temperature. The period of the compactified time, /?, is directly related to the 
temperature of the system in Lorentzian signature, through (feg = 1) 

T = j . (11.20) 

If we assume the standard Hawking temperature formula, T — kh/2tt, the period /3 in (|II.19|> can 
be re-written according to 

y/V'(f)W'{f) 

K H = ~ g ' (IL2i) 

which for the metric (|II.2[) reads kh — \al '(r) Hb(r) jj . Comparison between the latter and (|II.12|) 
fixes the normalization of the Killing vector £ a to be A = h. What is most important to us is 
that kh 7^ 0, so that we may conclude that the Killing horizon is of the bifurcate type. We may 
anticipate that this is not the unique surface gravity which can defined for a generic spherically 
symmetric static black hole. For the sake of simplicity, we shall postpone this discussion to the 
Conclusion an alternative definition. 

Given these preliminary remarks, we are now in the position to apply Wald's argument |30) to 
derive the black hole entropy associated to the Killing horizon of the solution (111.71) . 
Following [30l [31j [34], the explicit calculation of the black hole entropy Sw of the horizon 
r = rji = (c/p) p is provided by the formula 

ffJSf ^ (0) 



S w = ~2n<p — e ab e cd Jh {2) d9 dcp , (11.22) 

J *= ^H t \oR abcd J V 

where «Sf = J£(R a bcd, gab, V ' a Rbcde, ■ • ■ ) is the Lagrangian density of any general theory of gravity, 
in the specific case, 

^ (Rabcd, gab, V aRbcde,.-.) = (R + VSaVC?) . (11.23) 

107T 

The hatted variable, e a b, is the binormal vector to the (bifurcate) horizon: it is antisymmetyric 
under the exchange of a <H> b and normalized so that e a bi ab = —2. For the metric (|II.2[) . the 
binormal turns out to be 

e ab = b(r)(6° a 5t-5 l a 6 b ). (11.24) 

The induced volume form on the bifurcate surface r = ru,t ^constant is represented by \Jh(o\ dO d(f> 
where, for any spherically symmetric metric, \/h(2) — r 2 sm.9 and the angular variables 0, 4> run 
over the intervals [0,7r], [0, 2%), respectively. 

Finally, the superscript (0) indicates that the partial derivative 8J£ / 5R abc d is evaluated on shell. 
The variation of the Lagrangian density with respect to R abc d is performed as if R a bcd and the 
metric g ab are independent. 
In the specific case, equation pi.22[) becomes 

/ §<?> \ (°) 

S w = - W ff b 2 (r H ) — , (11.25) 



\5Rqioi . 

with &/h the area of the black hole horizon. Let us compute the Lagrangian variation, 
16tt (tfJSf) = SR + V3 a SiVc 2 ) 

= \{g ac g hd - g ad g bc )5Rabcd + ^{C 2 )--- S(C 2 ) . (11.26) 
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Using the fact that C 2 = R a bcdR abcd - 2R ab R ab + \R 2 , we get, 

S ^ —\-(n ac n bd a ad n bc \ + ^?-(r 2 \-h . 

167T 2 (ff 9 ~ 9 9 )+ 2 (C } 



SR 



abed 



2R abccl ~ {g ac R bd + g ba R ac - g ad R bc - g bc R ad ) + -(g ac g bd - g aa g bc )R 



11.27) 



In the specific, 



0101 



I 

327 



. 9 °y 1 + 



2jR 0101 _ g 0i? ll _ 9 H^00 + LgMgllR 



H 



Since in general, tr C n = (-§)" [2 + (-2) 2 -™]X™, for n > and 



X(r) = ^[r 2 a" + 2(a - 1) - 2ra'} + — [Zra'b' - 2a(b' - rb"" 



rb 



for the metric (|II.2[) . we may write 



r&\ H = -L \[r 2 a" + 2{a - 1) - 2m'] + \[3ra'b' - 2a(b' - rb")} 

rb 



V3 



H 



(11.28) 



(11.29) 



(11.30) 



Taking together (fTL25ll . (|IL28|| and (|IL30| . for both the solutions ((IT~5|) and (|ITT7|) . we finally have 
that the horizon entropy for the Deser et al. black hole is 



Sw = (1 + ea) 



where s := 



+1, 

-1, a>\ 



(11.31) 




fi 



According to (|II.31|) the entropy predicted by Wald's formula restricts considerably the space 
of the <7 parameter with respect to our previous considerations. In fact, as shown by Figure 1, the 
entropy of the black hole is positive only as far as <r G (— 1, 4]. For a = — 1, the entropy vanishes 
suggesting (but we leave this to the level of a speculation) that, for this value of a, the number 
of microscopic configurations realizing the black hole is only one. For a G (—1,0), the entropy of 
Deser's black hole is always smaller than its value in general relativity. Notice also that for a — 4, 
the entropy function is continuous even if the black hole metric changes. However, as pointed out 
above, such solution is not physical because of its pathological asymptotic behaviour. 
En passant, we notice how Wald's entropy could be computed equally well following [34J. Intro- 
ducing a new radial co-ordinate p such that 

p(r):= rv (11.32) 



the metric (III. 71) transforms to 



with 



ds 2 = -h{p)dt 2 + ^- + q(p)dfl 2 (11.33) 

Hp) 



Hp) = (^) 2(1 P) (p- —Lr J > ?G>) = (pk^p) 2p - (n.34) 

V pk " pj 



k 

This time, Wald's entropy (|II.31[) will follow from 



(0) 



S w = -torj> r= ^ [jji^J q(p)dfl 2 . (11.35) 

III The Clifton-Barrow solution 

The Clifton-Barrow solution starts from the following modified-gravity action (evaluated in the 
vacuum space): 

I CB = / d 4 x^ (—) ■ (HI.36) 
Jm \ X J 

Here, S is a constant and x is a dimensional parameter. We can choose \ = 167rG 1+<5 . When 
8 = 0, we recover the Hilbert-Einstein action of General Relativity. 

Taking the variation of the action with respect to the metric g^, we obtain: 

(d°d r R „ £ .d a Rd T R\ ( 1 + 2(5 \ , ttto . 

^ = \ — R ' — & — ) I 9 ^ a9uT + 2(l-8) 9 ^ 9aT ) ' ( m -37) 

Looking for static, spherically symmetric metric of the type, 

dr 2 

ds 2 = -V{r)dt 2 + -— + r 2 dft 2 , (111.38) 
W(r) 

we find the Clifton-Barrow solution of Equation (|III.37|) : 



T 



25(l+2(5)/(l-(5) 



V ^ =[- ) [ 1 + ' ( IIL39 ) 

W ^ = (l-25 + 46 2 )(l- 26 - 25 2 ) (* + r (i-2i+4^)/(i-«)) ' (IIL40) 



C and ro are dimensional constants. 
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In a similar way with respect to the previous section, we can see that the hypersurface 
r = r H := (_C)( 1 -< 5 )/( 1 - 2 < 5+4 < 5 ), for which W(r H ) = and d r W(r H ) ^ 0, determines an event 
horizon, and, since C < 0, the Clifton-Barrow metric is a Black Hole solution. 

According to Equation (|II.21|) . we recover the Killing-horizon surface gravity 



KH 




25 + 4<5 2 ) r£ 



(2<5+2<5 2 -l)/(l-<5) 



25 - 25 2 ) r Hws)/(i-S) 



(111.41) 



which can be used to find the Hawking temperature T — kh/2tt. 

As a last remark, we are able to derive the Black Hole entropy associated to the event horizon 
of the Clifton-Barrow solution. For modified gravity F(i?)-theories (where the gravity lagrangian 
is a function F(R) of the Ricci scalar only), it is easy to see that the Wald formula in Equation 
pi.22[) is simplified as 

dF(R) 



dR 



In our case, F(R) = R 1+s /x, so we find 



S W = ^(1 

X 



5) 



65(1 + 5) 



i <5 



(25 2 + 25 - l)r 2 H 



(111.42) 



(111.43) 



proved by the fact that on the Clifton-Barrow solution R = 65(1 + 5)/((25 2 + 25 — l)?" 2 ). 

In order to have the positive sign of entropy, we must require 5 > (V3 — l)/2. The solutions 
with < 5 < (y/3 — l)/2 are unphysical, whereas for 5 — we find the result of General Relativity. 



IV Conclusions 

Despite the great success enjoyed by modified theories of gravity, we have seen that only two 
non-trivial, static, spherically symmetric, vacuum black-hole solutions are known so far. Their 
thermodynamic properties have been taken into considerations. We have shown that the solutions 
we considered in §2 and 3 possess a Killing horizon with a Killing vector £ a ~ 9° associated which 
cannot be defined unambiguously due to the fact that the spacetimes are not asymptotically flat. 
What is most important, however, is that we have been able to deduce a non- vanishing Killing 
temperature for such horizons. Of course, this is not the only temperature we can define for such 
horizons. As shown in [36J, in spherically symmetric spacetimes always exists a Kodama vector 
field K whose defining property is that (G a bK b )' ,a = 0. The Kodama vector turns out to be 
time/light/space-like in untrapped/marginal/ trapped spacetime regions; it gives a preferred flow 
of time generalizing the Killing time flow familiar to static cases; it makes possible to define an 
invariant particle energy even in non-stationary spacetimes and it associates a Kodama-Hayward 
surface gravity to any future outer trapping horizon |37| . In static, asymptotically flat spacetimes, 
both the Killing and Kodama vectors coincide, so that they give rise to the same concepts of energy 
and temperature. In static, non-asymptotically flat spaces, they are both ambiguous and can differ 
by normalizations, but nonetheless the ratio "energy /surface gravity" remains fixed [B"51 14"0] , 
This means that as far as the Killing temperature associated with the black holes mentioned here 
is non- vanishing, also their Kodama-Hayward temperature will be so. On the other hand, that the 
horizons we are concerned are of the bifurcate type means we are in Wald's hypothesis in order 
to compute their entropy. In this sense, equations (jII.3lHTlI.43p and Figure 1 represent our main 
results. 

To complete the picture of thermodynamic features of black holes in modified theories of 
gravity, it would be necessary to formulate a consistent definition of their mass. As it is well 
known, in modified theories of gravity the first law of thermodynamics generally requires a work 
term even in vacuum solutions something which makes the first law useless in the situations at 
hand. Quite recently some attempts have been put forward in order to answer the question, but 
only for asymptotically flat spacetimes, cf. [4"Ti [4"2] . 
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In principle, a powerful tool to evaluate the black hole mass in a theory of the type Jz? = R + (. . . ) 
is represented by the so-called Brown- York quasi- local mass [15114^1143] . In static, spherically sym- 
metric spacetimes where the metric can be put in the form (|II.2[) the BY mass reads 



Mby = ra(r)b(r) 



z(°)(r) 
a(r) 



(IV.44) 



with a'°)(r) an arbitrary function which determines the zero of the energy for a background 
spacetime and r is the radius of the spacelike hypersurface boundary. When the spacetime is 
asymptotically flat, the ADM mass M is the Mby determined in pV.44|) in the limit r — > oo. 
If no cosmological horizon is present, the large r limit of (|IV.44p is used to determine the mass. 
However, this approach is known to fail whenever the matter action (i.e. what we have represented 
with (...) few line above) contains derivatives of the metric as it is the case of the Deser et al. 



action, (|II.l[l . 

Another quasi-local energy definition well known in general relativity and fully employed in 
spherical symmetry is the so-called Misner-Sharp energy |46| which can be proved to be the 
conserved charge generated by the Kodama vector K [4"Tl 148|. In the last few years, different 
authors have tried to generalize the Misner-Sharp energy definition to wider classes of gravity 
theories [49 , 50J . But even if Cai et al. provide a general formula for the generalized MS energy in 
f(R) gravity, this does not produce any explicit, useful, result for the Clifton-Barrow black hole. 

In conclusion, we have succeeded in computing two of three most relevant thermodynamic 
parameters (entropy and temperature) of the known black hole solutions in modified theories of 
gravity; the mass resisting up to now to any attack led by conventional methods. 
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